The genome of the rayed Mediterranean limpet Patella caerulea (Linnaeus, 1758).

Patella caerulea (Linnaeus, 1758) is a mollusc limpet species of the class Gastropoda. Endemic to the Mediterranean Sea, it is considered a keystone species due to its primary role in structuring and regulating the ecological balance of tidal and subtidal habitats. It is currently being used as a bioindicator to assess the environmental quality of coastal marine waters and as a model species to understand adaptation to ocean acidification. Here, we provide a high-quality reference genome assembly and annotation for P. caerulea. We generated ∼30 Gb of Pacific Biosciences high-fidelity data from a single individual and provide a final 749.8 Mb assembly containing 62 contigs, including the mitochondrial genome (14,938 bp). With an N50 of 48.8 Mb and 98% of the assembly contained in the 18 largest contigs, this assembly is near chromosome-scale. Benchmarking Universal Single-Copy Orthologs scores were high (Mollusca, 87.8% complete; Metazoa, 97.2% complete) and similar to metrics observed for other chromosome-level Patella genomes, highlighting a possible bias in the Mollusca database for Patellids. We generated transcriptomic Illumina data from a second individual collected at the same locality and used it together with protein evidence to annotate the genome. A total of 23,938 protein-coding gene models were found. By comparing this annotation with other published Patella annotations, we found that the distribution and median values of exon and gene lengths was comparable with other Patella species despite different annotation approaches. The present high-quality P. caerulea reference genome, available on GenBank (BioProject: PRJNA1045377; assembly: GCA_036850965.1), is an important resource for future ecological and evolutionary studies.

Single-cell transcriptomics of the human parasite Schistosoma mansoni first intra-molluscan stage reveals tentative tegumental and stem-cell regulators.

Schistosomiasis is a major Neglected Tropical Disease, caused by the infection with blood flukes in the genus Schistosoma. To complete the life cycle, the parasite undergoes asexual and sexual reproduction within an intermediate snail host and a definitive mammalian host, respectively. The intra-molluscan phase provides a critical amplification step that ensures a successful transmission. However, the cellular and molecular mechanisms underlying the development of the intra-molluscan stages remain poorly understood. Here, single cell suspensions from S. mansoni mother sporocysts were produced and sequenced using the droplet-based 10X Genomics Chromium platform. Six cell clusters comprising two tegument, muscle, neuron, parenchyma and stem/germinal cell clusters were identified and validated by in situ hybridisation. Gene Ontology term analysis predicted key biological processes for each of the clusters, including three stem/germinal sub-clusters. Furthermore, putative transcription factors predicted for stem/germinal and tegument clusters may play key roles during parasite development and interaction with the intermediate host.

Revisiting the evolution of Family B1 GPCRs and ligands: Insights from mollusca.

Family B1 G protein-coupled receptors (GPCRs) are one of the most well studied neuropeptide receptor families since they play a central role in many biological processes including endocrine, gastrointestinal, cardiovascular and reproduction in animals. The genes for these receptors emerged from a common ancestral gene in bilaterian genomes and evolved via gene/genome duplications and deletions in vertebrate and invertebrate genomes. Their existence and function have mostly been characterized in vertebrates and few studies exist in invertebrate species. Recently, an increased interest in molluscs, means a series of genomes have become available, and since they are less modified than insect and nematode genomes, they are ideal to explore the origin and evolution of neuropeptide gene families. This review provides an overview of Family B1 GPCRs and their peptide ligands and incorporates new data obtained from Mollusca genomes and taking a comparative approach challenges existing models on their origin and evolution.

Phylogenomic reconstruction of Solenogastres (Mollusca, Aplacophora) informs hypotheses on body size evolution.

Body size is a fundamental characteristic of animals that impacts every aspect of their biology from anatomical complexity to ecology. In Mollusca, Solenogastres has been considered important to understanding the group's early evolution as most morphology-based phylogenetic reconstructions placed it as an early branching molluscan lineage. Under this scenario, molluscs were thought to have evolved from a small, turbellarian-like ancestor and small (i.e., macrofaunal) body size was inferred to be plesiomorphic for Solenogastres. More recently, phylogenomic studies have shown that aplacophorans (Solenogastres + Caudofoveata) form a clade with chitons (Polyplacophora), which is sister to all other molluscs, suggesting a relatively large-bodied (i.e., megafaunal) ancestor for Mollusca. Meanwhile, recent investigations into aplacophoran phylogeny have called the assumption that the last common ancestor of Solenogastres was small-bodied into question, but sampling of meiofaunal species was limited, biasing these studies towards large-bodied taxa and leaving fundamental questions about solenogaster body size evolution unanswered. Here, we supplemented available data with transcriptomes from eight diverse meiofaunal species of Solenogastres and conducted phylogenomic analyses on datasets of up to 949 genes. Maximum likelihood analyses support the meiofaunal family Meiomeniidae as the sister group to all other solenogasters, congruent with earlier ideas of a small-bodied ancestor of Solenogastres. In contrast, Bayesian Inference analyses support the large-bodied family Amphimeniidae as the sister group to all other solenogasters. Investigation of phylogenetic signal by comparing site-wise likelihood scores for the two competing hypotheses support the Meiomeniidae-first topology. In light of these results, we performed ancestral character state reconstruction to explore the implications of both hypotheses on understanding of Solenogaster evolution and review previous hypotheses about body size evolution and its potential consequences for solenogaster biology. Both hypotheses imply that body size evolution has been highly dynamic over the course of solenogaster evolution and that their relatively static body plan has successfully allowed for evolutionary transitions between meio-, macro- and megafaunal size ranges.

A chromosome-level genome for the nudibranch gastropod Berghia stephanieae helps parse clade-specific gene expression in novel and conserved phenotypes.

BACKGROUND: How novel phenotypes originate from conserved genes, processes, and tissues remains a major question in biology. Research that sets out to answer this question often focuses on the conserved genes and processes involved, an approach that explicitly excludes the impact of genetic elements that may be classified as clade-specific, even though many of these genes are known to be important for many novel, or clade-restricted, phenotypes. This is especially true for understudied phyla such as mollusks, where limited genomic and functional biology resources for members of this phylum have long hindered assessments of genetic homology and function. To address this gap, we constructed a chromosome-level genome for the gastropod Berghia stephanieae (Valdés, 2005) to investigate the expression of clade-specific genes across both novel and conserved tissue types in this species. RESULTS: The final assembled and filtered Berghia genome is comparable to other high-quality mollusk genomes in terms of size (1.05 Gb) and number of predicted genes (24,960 genes) and is highly contiguous. The proportion of upregulated, clade-specific genes varied across tissues, but with no clear trend between the proportion of clade-specific genes and the novelty of the tissue. However, more complex tissue like the brain had the highest total number of upregulated, clade-specific genes, though the ratio of upregulated clade-specific genes to the total number of upregulated genes was low. CONCLUSIONS: Our results, when combined with previous research on the impact of novel genes on phenotypic evolution, highlight the fact that the complexity of the novel tissue or behavior, the type of novelty, and the developmental timing of evolutionary modifications will all influence how novel and conserved genes interact to generate diversity.

A chromosome-level genome assembly of a deep-sea symbiotic Aplacophora mollusc Chaetoderma sp.

The worm-shaped, shell-less Caudofoveata is one of the least known groups of molluscs. As early-branching molluscs, the lack of high-quality genomes hinders our understanding of their evolution and ecology. Here, we report a high-quality chromosome-scale genome of Chaetoderma sp. combining PacBio, Illumina, and high-resolution chromosome conformation capture sequencing. The final assembly has a size of 2.45 Gb, with a scaffold N50 length of 141.46 Mb, and is anchored to 17 chromosomes. Gene annotations showed a high level of accuracy and completeness, with 23,675 predicted protein-coding genes and 94.44% of the metazoan conserved genes by BUSCO assessment. We further present 16S rRNA gene amplicon sequencing of the gut microbiota in Chaetoderma sp., which was dominated by the chemoautotrophic bacteria (phylum Gammaproteobacteria). This chromosome-level genome assembly presents the first genome for the Caudofoveata, which constitutes an important resource for studies ranging from molluscan evolution, symposium, to deep-sea adaptation.

A combined phylogenetic strategy illuminates the evolution of Goniodorididae nudibranchs (Mollusca, Gastropoda, Heterobranchia).

Goniodorididae is a family of small dorid nudibranchs distributed worldwide that feed on entoprocts, ascidians, and bryozoans. The evolutionary relationships between its taxa have been uncertain due to the limited taxa available for phylogenetic analyses; some genera being paraphyletic. The family includes a remarkable number of synonymized genera in which the species richness is unequally distributed, while some genera have dozens of species others are monospecific. Some clades are very uniform morphologically while others are considered highly variable. To increase backbone phylogenetic resolution a target enrichment approach of ultra-conserved elements was aimed at representative Goniodorididae species for the first time. Additionally, we increase species representation by including mitochondrial markers cytochrome c oxidase subunit I and ribosomal RNA 16S as well as nuclear Histone 3 and ribosomal RNA 18S from 109 Goniodorididae species, out of approximately 160 currently valid species. Maximum likelihood and Bayesian inference analyses were performed to infer the phylogeny of the family. As a result, two subfamilies and eleven genera were elucidated. The synonymized genera Bermudella, Cargoa, and Ceratodoris are here resurrected and a new genus, Naisdoris gen. nov., is described. The clades included taxa with shared prey preference, showing that trophic behavior could have driven species evolution and morphological uniqueness within the family Goniodorididae.

Molluscan Genomes Reveal Extensive Differences in Photopigment Evolution Across the Phylum.

In animals, opsins and cryptochromes are major protein families that transduce light signals when bound to light-absorbing chromophores. Opsins are involved in various light-dependent processes, like vision, and have been co-opted for light-independent sensory modalities. Cryptochromes are important photoreceptors in animals, generally regulating circadian rhythm, they belong to a larger protein family with photolyases, which repair UV-induced DNA damage. Mollusks are great animals to explore questions about light sensing as eyes have evolved multiple times across, and within, taxonomic classes. We used molluscan genome assemblies from 80 species to predict protein sequences and examine gene family evolution using phylogenetic approaches. We found extensive opsin family expansion and contraction, particularly in bivalve xenopsins and gastropod Go-opsins, while other opsins, like retinochrome, rarely duplicate. Bivalve and gastropod lineages exhibit fluctuations in opsin repertoire, with cephalopods having the fewest number of opsins and loss of at least 2 major opsin types. Interestingly, opsin expansions are not limited to eyed species, and the highest opsin content was seen in eyeless bivalves. The dynamic nature of opsin evolution is quite contrary to the general lack of diversification in mollusk cryptochromes, though some taxa, including cephalopods and terrestrial gastropods, have reduced repertoires of both protein families. We also found complete loss of opsins and cryptochromes in multiple, but not all, deep-sea species. These results help set the stage for connecting genomic changes, including opsin family expansion and contraction, with differences in environmental, and biological features across Mollusca.

Genomic Insights into Mollusk Terrestrialization: Parallel and Convergent Gene Family Expansions as Key Facilitators in Out-of-the-Sea Transitions.

Animals abandoned their marine niche and successfully adapted to life on land multiple times throughout evolution, providing a rare opportunity to study the mechanisms driving large scale macroevolutionary convergence. However, the genomic factors underlying this process remain largely unknown. Here, we investigate the macroevolutionary dynamics of gene repertoire evolution during repeated transitions out of the sea in mollusks, a lineage that has transitioned to freshwater and terrestrial environments multiple independent times. Through phylogenomics and phylogenetic comparative methods, we examine ∼100 genomic data sets encompassing all major molluskan lineages. We introduce a conceptual framework for identifying and analyzing parallel and convergent evolution at the orthogroup level (groups of genes derived from a single ancestral gene in the species in question) and explore the extent of these mechanisms. Despite deep temporal divergences, we found that parallel expansions of ancient gene families played a major role in facilitating adaptation to nonmarine habitats, highlighting the relevance of the preexisting genomic toolkit in facilitating adaptation to new environments. The expanded functions primarily involve metabolic, osmoregulatory, and defense-related systems. We further found functionally convergent lineage-exclusive gene gains, while family contractions appear to be driven by neutral processes. Also, genomic innovations likely contributed to fuel independent habitat transitions. Overall, our study reveals that various mechanisms of gene repertoire evolution-parallelism, convergence, and innovation-can simultaneously contribute to major evolutionary transitions. Our results provide a genome-wide gene repertoire atlas of molluskan terrestrialization that paves the way toward further understanding the functional and evolutionary bases of this process.

Hiker, a new family of DNA transposons encoding transposases with DD35E motifs, displays a distinct phylogenetic relationship with most known DNA transposon families of IS630-Tc1-mariner (ITm).

DNA transposons play a crucial role in determining the size and structure of eukaryotic genomes. In this study, a new family of IS630-Tc1-mariner (ITm) DNA transposons, named Hiker (HK), was identified. HK is characterized by a DD35E catalytic domain and is distinct from all previously known families of the ITm group. Phylogenetic analyses showed that DD35E/Hiker forms a monophyletic clade with DD34E/Gambol, indicating that they may represent a separate superfamily of ITm. A total of 178 Hiker species were identified, with 170 found mainly in Actinopterygii, one in Chondrichthyes, six in Anura and one in Mollusca. Gambol (GM), on the other hand, are found in invertebrates, with 18 in Arthropoda and one in Platyhelminthes. Hiker transposons have a total length ranging from 2.14 to 3.67 kb and contain a single open reading frame that encodes a protein of approximately 370 amino acids (range 311-413 aa). They are flanked by short terminal inverted repeats (TIRs) of 16-30 base pairs and two base pair (TA) target-site duplications. In contrast, most transposons of the Gambol family have a total length of 1.35-5.96 kb, encode a transposase protein of approximately 350 amino acids (range 306-374 aa), and are flanked by TIRs that range from 32 to 1097 bp in length. Both Hiker and Gambol transposases have several conserved motifs, including helix-turn-helix (HTH) motifs and a DDE domain. Our study observed multiple amplification waves and repeated horizontal transfer (HT) events of HK transposons in vertebrate genomes, indicating their role in diversifying and shaping the genomes of Actinopterygii, Chondrichthyes, and Anura. Conversely, GM transposons showed few Horizontal transfer events. According to cell-based transposition assays, most HK transposons are likely inactive due to the truncated DNA binding domains of their transposases. We present an updated classification of the ITm group based on these findings, which will enhance the understanding of both the evolution of ITm transposons and that of their hosts.

Role of maternal spiralian-specific homeobox gene SPILE-E in the specification of blastomeres along the animal-vegetal axis during the early cleavage stages of mollusks.

Spiralians, one of the major clades of bilaterians, share a unique development known as spiralian development, characterized by the formation of tiers of cells called quartets, which exhibit different developmental potentials along the animal-vegetal axis. Recently, spiralian-specific TALE-type homeobox genes (SPILE) have been identified, some of which show zygotic and staggered expression patterns along the animal-vegetal axis and function in quartet specification in mollusks. However, it is unclear which maternal molecular components control the zygotic expression of these transcription factors. In this study, we focused on SPILE-E, a maternal transcription factor, and investigated its expression and function in mollusks. We found that the maternal and ubiquitous expression of SPILE-E in the cleavage stages is conserved in molluskan species, including limpets, mussels, and chitons. We knocked down SPILE-E in limpets and revealed that the expression of transcription factors specifically expressed in the first quartet (1q2 ; foxj1b) and second quartet (2q; SPILE-B) was abolished, whereas the macromere-quartet marker (SPILE-C) was ectopically expressed in 1q2 in SPILE-E morphants. Moreover, we showed that the expression of SPILE-A, which upregulates SPILE-B but represses SPILE-C expression, decreased in SPILE-E morphants. Consistent with changes in the expression pattern of the above transcription factors, SPILE-E-morphant larvae exhibited patchy or complete loss of expression of marker genes of ciliated cells and shell fields, possibly reflecting incomplete specification of 1q2 and 2q. Our results provide a molecular framework for quartet specification and highlight the importance of maternal lineage-specific transcription factors in the development and evolution of spiralians.

Mosquito (MS), a DD37E Family of Tc1/Mariner, Displaying a Distinct Evolution Profile from DD37E/TRT and DD37E/L18.

Diverse Tc1/mariner elements with the DD37E signature have been detected. However, their evolutionary relationship and profiles are largely unknown. Using bioinformatics methods, we defined the evolution profile of a Tc1/Mariner family, which harbors the catalytic domain with the DD37E signature, and renamed it DD37E/Mosquito (MS). MS transposons form a separate monophyletic clade in the phylogenetic tree, distinct from the other two groups of elements with the DD37E signature, DD37E/L18 and DD37E/TRT (transposon related to Tc1), and represent a very different taxonomic distribution from that of DD37E/TRT. MS is only detected in invertebrate and is mostly present in Arthropoda, as well as in Cnidaria, Ctenophora, Mollusca, Nematoda, and Platyhelminthes, with a total length of about 1.3 kb, containing an open reading frame (ORF) encoding about 340 amino acids transposases, with a conserved DD37E catalytic domain. The terminal inverted repeat (TIR) lengths range from 19 bp to 203 bp, and the target site duplication (TSD) is TA. We also identified few occurrences of MS horizontal transfers (HT) across lineages of diptera. In this paper, the distribution characteristics, structural characteristics, phylogenetic evolution, and horizontal transfer of the MS family are fully analyzed, which is conducive to supplementing and improving the Tc1/Mariner superfamily and excavating active transposons.

Proton channels in molluscs: A new bivalvian-specific minimal HV 4 channel.

Recently, three proton channels (HV ) have been identified and characterized in Aplysia californica (AcHV 1-3). Focusing on AcHV 1 and AcHV 2, analysis of Transcriptome Shotgun Assembly and genomic databases of 91 molluscs identified HV homologous channels in other molluscs: channels homologous to AcHV 1 and to AcHV 2 were found in 90 species (56 full-length sequences) and in 33 species (18 full-length sequences), respectively. Here, we report the discovery of a fourth distinct proton channel family, HV 4. This new family has high homology to AcHV 1 and AcHV 2 and was identified only in bivalvian molluscs (13 species, 12 full-length sequences). Typically, these channels possess an extracellular S1-S2 loop of intermediate size (~ 20 amino acids) compared to the shorter loops of molluscan HV 1 channels (~ 13 amino acids) and the much larger loops of molluscan HV 2 channels (> 65 amino acids). The characteristic voltage-sensor motif in S4 possesses only two arginine residues with the common third arginine being replaced by a lysine. Moreover, HV 4 channels are much smaller with only around 200 amino acids in total length. The smallest functional channel found so far in nature (189 amino acids) is expressed in the pacific oyster Crassostrea gigas (CgHV 4) and might be considered an archetypical minimal proton channel. Functional expression and electrophysiological characterization demonstrated that CgHV 4 shares distinctive hallmarks of other investigated proton channels as high proton selectivity, slow activation, and pH- and voltage-regulated gating. This work is the first description of a HV 4 type channel, adding a new member to the recently expanded family of proton channels.

The Complete Mitochondrial Genome and Gene Arrangement of the Enigmatic Scaphopod Pictodentalium vernedei.

The enigmatic scaphopods, or tusk shells, are a small and rare group of molluscs whose phylogenomic position among the Conchifera is undetermined, and the taxonomy within this class also needs revision. Such work is hindered by there only being a very few mitochondrial genomes in this group that are currently available. Here, we present the assembly and annotation of the complete mitochondrial genome from Dentaliida Pictodentalium vernedei, whose mitochondrial genome is 14,519 bp in size, containing 13 protein-coding genes, 22 tRNA genes and two rRNA genes. The nucleotide composition was skewed toward A-T, with a 71.91% proportion of AT content. Due to the mitogenome-based phylogenetic analysis, we defined P. vernedei as a sister to Graptacme eborea in Dentaliida. Although a few re-arrangements occurred, the mitochondrial gene order showed deep conservation within Dentaliida. Yet, such a gene order in Dentaliida largely diverges from Gadilida and other molluscan classes, suggesting that scaphopods have the highest degree of mitogenome arrangement compared to other molluscs.

Application of palaeogenetic techniques to historic mollusc shells reveals phylogeographic structure in a New Zealand abalone.

Natural history collections worldwide contain a plethora of mollusc shells. Recent studies have detailed the sequencing of DNA extracted from shells up to thousands of years old and from various taphonomic and preservational contexts. However, previous approaches have largely addressed methodological rather than evolutionary research questions. Here, we report the generation of DNA sequence data from mollusc shells using such techniques, applied to Haliotis virginea Gmelin, 1791, a New Zealand abalone, in which morphological variation has led to the recognition of several forms and subspecies. We successfully recovered near-complete mitogenomes from 22 specimens including 12 dry-preserved shells up to 60 years old. We used a combination of palaeogenetic techniques that have not previously been applied to shell, including DNA extraction optimized for ultra-short fragments and hybridization-capture of single-stranded DNA libraries. Phylogenetic analyses revealed three major, well-supported clades comprising samples from: (1) The Three Kings Islands; (2) the Auckland, Chatham and Antipodes Islands; and (3) mainland New Zealand and Campbell Island. This phylogeographic structure does not correspond to the currently recognized forms. Critically, our nonreliance on freshly collected or ethanol-preserved samples enabled inclusion of topotypes of all recognized subspecies as well as additional difficult-to-sample populations. Broader application of these comparatively cost-effective and reliable methods to modern, historical, archaeological and palaeontological shell samples has the potential to revolutionize invertebrate genetic research.

Hemoglobin wonders: a fascinating gas transporter dive into molluscs.

Hemoglobin (Hb) has been identified in at least 14 molluscan taxa so far. Research spanning over 130 years on molluscan Hbs focuses on their genes, protein structures, functions, and evolution. Molluscan Hbs are categorized into single-, two-, and multiple-domain chains, including red blood cell, gill, and extracellular Hbs, based on the number of globin domains and their respective locations. These Hbs exhibit variation in assembly, ranging from monomeric and dimeric to higher-order multimeric forms. Typically, molluscan Hbs display moderately high oxygen affinity, weak cooperativity, and varying pH sensitivity. Hb's potential role in antimicrobial pathways could augment the immune defense of bivalves, which may be a complement to their lack of adaptive immunity. The role of Hb as a respiratory protein in bivalves likely originated from the substitution of hemocyanin. Molluscan Hbs demonstrate adaptive evolution in response to environmental changes via various strategies (e.g. increasing Hb types, multimerization, and amino acid residue substitutions at key sites), enhancing or altering functional properties for habitat adaptation. Concurrently, an increase in Hb assembly diversity, coupled with a downward trend in oxygen affinity, is observed during molluscan differentiation and evolution. Hb in Protobranchia, Heteroconchia, and Pteriomorphia bivalves originated from separate ancestors, with Protobranchia inheriting a relative ancient molluscan Hb gene. In bivalves, extracellular Hbs share a common origin, while gill Hbs likely emerged from convergent evolution. In summary, research on molluscan Hbs offers valuable insights into the origins, biological variations, and adaptive evolution of animal Hbs.

Molecular Phylogenetics of Physa acuta (Pulmonata: Basommatophora): an invasive species in Central Punjab Pakistan / Molecular Phylogenetics of Physa acuta (Pulmonata: Basommatophora): uma espécie invasora no centro de Punjab, Paquistão

Physids belong to Class Gastropoda; belong to Phylum Mollusca and being bioindicators, intermediate hosts of parasites and pests hold a key position in the ecosystem. There are three species of Genus Physa i.e. P. fontinalis, Physa acuta and P. gyrina water bodies of Central Punjab and were characterized on the basis of molecular markers High level of genetic diversity was revealed by polymorphic RAPD, however SSR markers were not amplified. The multivariate analysis revealed polymorphism ranging from 9.09 percent to 50 percent among the three Physid species. Total number of 79 loci were observed for the three species under study and 24 loci were observed to be polymorphic. These RAPD fragment(s) can be developed into co dominant markers (SCAR) by cloning and can be further sequenced for the development of the Physa species specific markers to identify the introduced and native species in Pakistan.

Os físidos pertencem à classe Gastropoda; pertencem ao filo Mollusca e, sendo bioindicadores, hospedeiros intermediários de parasitas e pragas, ocupam uma posição-chave no ecossistema. Existem três espécies do gênero Physa, ou seja, P. fontinalis, Physa acuta e P. gyrina em corpos d’água do Punjab Central e foram caracterizadas com base em marcadores moleculares. Alto nível de diversidade genética foi revelado por RAPD polimórfico, no entanto os marcadores SSR não foram amplificados. A análise multivariada revelou polimorfismo variando de 9,09% a 50% entre as três espécies de Physid. Um número total de 79 loci foi observado para as três espécies em estudo e 24 loci foram observados como polimórficos. Esses fragmentos RAPD podem ser desenvolvidos em marcadores codominantes (SCAR) por clonagem e podem ser posteriormente sequenciados para o desenvolvimento de marcadores específicos da espécie Physa para identificar as espécies introduzidas e nativas no Paquistão.

Exploring the Potential of Metatranscriptomics to Describe Microbial Communities and Their Effects in Molluscs.

Metatranscriptomics has emerged as a very useful technology for the study of microbiomes from RNA-seq reads. This method provides additional information compared to the sequencing of ribosomal genes because the gene expression can also be analysed. In this work, we used the metatranscriptomic approach to study the whole microbiome of mussels, including bacteria, viruses, fungi, and protozoans, by mapping the RNA-seq reads to custom assembly databases (including the genomes of microorganisms publicly available). This strategy allowed us not only to describe the diversity of microorganisms but also to relate the host transcriptome and microbiome, finding the genes more affected by the pathogen load. Although some bacteria abundant in the metatranscriptomic analysis were undetectable by 16S rRNA sequencing, a common core of the taxa was detected by both methodologies (62% of the metatranscriptomic detections were also identified by 16S rRNA sequencing, the Oceanospirillales, Flavobacteriales and Vibrionales orders being the most relevant). However, the differences in the microbiome composition were observed among different tissues of Mytilus galloprovincialis, with the fungal kingdom being especially diverse, or among molluscan species. These results confirm the potential of a meta-analysis of transcriptome data to obtain new information on the molluscs' microbiome.

Molecular phylogeny of selected dorid nudibranchs based on complete mitochondrial genome.

Dorid nudibranchs are a large group of mollusks with approximately 2,000 recorded species. Although agreement exists on the monophyletic nature of the dorid nudibranch group, the interfamily relationships of the suborder are subject to debate. Despite efforts to elucidate this issue using short molecular markers, the conclusiveness of the findings has been hindered by branching polytomy. Mitogenomes are known to be effective markers for use in phylogenetic investigations. In this study, eight mitogenomes of dorid nudibranchs were decoded and analyzed. Gene content and structure showed little change among species, reflecting the conserved mitogenomes of dorid nudibranchs. For most genes, the direction was typical for nudibranchs; nevertheless, tRNACys had an inverse direction in Cadlinidae species. Phylogenetic trees based on nucleotide and amino acid datasets revealed a relatively consistent pattern of interfamily relationships with little difference for positions of Phyllidiidae and Cadlinidae. Species of Cadlinidae were clustered together and did not form a clade with Chromododidae. Additionally, Goniodorididae was sister to Aegiridae, whereas Discodoridae was sister to Dorididae. This finding was supported by tree topology test based on mitogenome data. The results of the present study indicate that complete mitogenomes are promising markers for investigating interfamily relationships among dorid nudibranchs.

Structural and evolutionary insights into the multidomain galectin from the red abalone Haliotis rufescens with specificity for sulfated glycans.

Galectins are an evolutionarily ancient family of lectins characterized by their affinity for ß-galactosides and a conserved binding site in the carbohydrate recognition domain (CRD). These lectins are involved in multiple physiological functions, including the recognition of glycans on the surface of viruses and bacteria. This feature supports their role in innate immune responses in marine mollusks. Here, we identified and characterized a galectin, from the mollusk Haliotis rufescens (named HrGal), with four CRDs that belong to the tandem-repeat type. HrGal was purified by affinity chromatography in a galactose-agarose resin and exhibited a molecular mass of 64.11 kDa determined by MALDI-TOF mass spectrometry. The identity of HrGal was verified by sequencing, confirming that it is a 555 amino acid protein with a mass of 63.86 kDa. This protein corresponds to a galectin reported in GenBank with accession number AHX26603. HrGal is stable in the presence of urea, reducing agents, and ions such as Cu2+ and Zn2+. The recombinant galectin (rHrGal) was purified from inclusion bodies in the presence of these ions. A theoretical model obtained with the AlphaFold server exhibits four non-identical CRDs, with a ß sandwich folding and the representative motifs for binding ß-galactosides. This allows us to classify HrGal within the tandem repeat galectin family. On the basis of a phylogenetic analysis, we found that the mollusk sequences form a monophyletic group of tetradomain galectins unrelated to vertebrate galectins. HrGal showed specificity for galactosides and glucosides but only the sulfated sugars heparin and ι-carrageenan inhibited its hemagglutinating activity with a minimum inhibitory concentration of 4 mM and 6.25 X 10-5% respectively. The position of the sulfate groups seemed crucial for binding, both by carrageenans and heparin.